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Abstract 

We study the detectabihty of the gravitational waves (GWs) from the Q-bah formation associated 
with the Affleck-Dine (AD) mechanism, taking into account both the dilution effects due to Q-ball 
domination and to finite temperature. The AD mechanism predicts the formation of nontopological 
solitons, Q-balls, from which GWs are generated. Q-balls with large conserved charge Q can 
produce a large amount of GWs. On the other hand, the decay rate of such Q-balls is so small that 
they may dominate the energy density of the Universe, which implies that GWs are significantly 
diluted and that their frequencies are redshifted during the Q-ball dominated era. Thus, the 
detectabihty of the GWs associated with the formation of Q-balls is determined by these two 
competing effects. We find that there is a finite but small parameter region where such GWs may 
be detected by future detectors such as DECIGO or BBO, only in the case when the thermal 
logarithmic potential dominates the potential of the AD field. Otherwise GWs from Q-balls would 
not be detectable even by these futuristic detectors: ^^q\y < 10"^^. Unfortunately, for such 
parameter region the present baryon asymmetry of the Universe can hardly be explained unless 
one fine-tunes A-terms in the potential. However the detection of such a GW background may give 
us an information about the early Universe, for example, it may suggest that the flat directions 
with B — L = are favored. 

PACS numbers: 98.80.Cq ; 04.30.Db 
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I. INTRODUCTION 



Primordial gravitational waves (GWs) provide us with a lot of important information 
about the early Universe because the interaction of GWs with matter is very weak and 
they carry the memory of cosmic history during and after inflation One of the main 
sources of such primordial GWs is cosmic inflation j2|, during which stochastic GWs are 
generated with a nearly scale-invariant spectrum with frequencies ranging from 10~^^ Hz 
to 10^ Hz, depending on inflation models. Though the amplitude is typically very small, 
next-generation gravitational detectors like DECIGO jsj and BBO j3] as well as indirect 
observations through the 5-modes in the cosmic microwave background (CMB) anisotropy 
[sl might be able to detect such GWs. The detection of such stochastic GWs can reveal the 
history of the early Universe such as the changes of the number of the massless degree of 
freedom [6|, the reheating temperature the neutrino decoupling [sj, the lepton asymmetry 
[qI, and so on. Additional GWs can also be generated from cosmological mechanisms like 
(local/global) phase transitions [lo[ and preheating after inflation 11, l2] or astrophysical 
origins like the collapse of Population HI stars fl3| . 

Recently, yet another interesting mechanism to produce GWs was proposed in Ref. [li 
in which it is shown that significant GWs are emitted during the Q-ball formation associ- 
ated with the Affleck-Dine (AD) mechanism of baryogenesis. In this scenario, the angular 
momentum, the rotation around the origin of scalar fields that carry the baryon or lepton 
number, is dynamically generated, which implies that the baryon or lepton asymmetry is 
produced. When initially almost homogeneous scalar fields start to rotate, they suffer from 
spatial instabilities, and their fluctuations begin to grow if the scalar field potential driving 
the field rotation is fiatter than the quadratic potential. Such fiuctuations finally settle down 
into nontopological solitons called Q-balls, whose existence and stability are guaranteed by 
a conserved charge, Q, associated with a global symmetry 15|. Since the formation pro- 



cess of such Q-balls is inhomogeneous and not spherical, GWs can be generated during the 
formation. 

In order to calculate the present amplitude of such GWs, one has to estimate not only the 
amount of GWs at the formation of Q-balls but also the dilution factor during the cosmic 



history after the production of the GWs. As shown in Ref. [IJ], the energy density of the 
GWs at the Q-ball formation is proportional to some powers of the field value of the Affleck- 
Dine condensate, which implies that the initial energy density of the GWs becomes large if 
the typical charge Q of Q-balls is large. On the other hand, the lifetime of Q-balls becomes 
longer for larger Q because the temperature at the decay of Q-balls is typically proportional 
to the inverse square root of the charge Q jl6|. Therefore, Q-balls with large Q can quickly 
dominate the energy density of the Universe and hence dilute the GWs significantly. Thus, 
the maximum value of the present amplitude of such GWs is determined by the balance of 
the above two competing effects. In Refs. [l^, in order to avoid the large dilution by the 
Q-ball dominance, Q-balls are assumed to decay quickly by some artificial effects, which are 
unnatural in the context of the minimal supersymmetric standard model (MSSM). Then, 
in this paper, we reconsider the decay of Q-balls without such effects and calculate the 
amplitude of GWs from Q-balls taking into account of the dilution factor correctly, which 
results in dramatic changes in the present amplitude of the GWs from the Q-ball formation. 

The properties of the Affleck-Dine mechanism [13] and the subsequent formation of Q- 
balls depend on the supersymmetry (SUSY) breaking mechanism 18, 19||. This is mainly 
because the effective potential of a flat direction and the gravitino mass are quite different for 
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different mediation mechanisms. There are many types of Q-balls such as gauge-mediation 



type [18|, |20|, gravity-mediation type [19|, |2l|, new type [22|, delayed type [23|, and so on. 



Moreover, the effective potential of a flat direction consists of thermal terms as well as the 
zero-temperature terms. The properties of Q-balls are quite different depending on which 
term dominates the potential energy. When thermal (logarithmic) effects |24l | dominate 



the effective potential, the formed Q-balls, called thermal log type |23|, have an interesting 



property. Namely, while the energy density of other types of Q-balls decreases like matter, 
that of the thermal log type Q-ball decreases at least as rapid as radiation. This is because 
the thermal logarithmic potential itself also decreases with the cosmic expansion while the 
number of Q-balls in a comoving volume does not change. Hence, this type of Q-balls cannot 
dominate the energy density of the Universe and do not dilute GWs. This is favorable for 
the detection of the GWs from the Q-ball formation because the dilution during the Q-ball 
dominated era is the main obstacle for the detection. In fact, as is explicitly shown in 
Appendix [Bl GWs from the Q-ball formation in the zero-temperature potential may not be 
detectable even by the next-generation detectors. Therefore, in this paper, we concentrate 
on the thermal log type Q-balls and estimate the present amplitudes and frequencies of the 
GWs at the formation of such Q-balls. We show that such GWs may be detected by the 
next-generation gravitational detectors like DECIGO and BBO if particular conditions of 
reheating temperature, the initial field value of the AD field, gravitino mass and messenger 
mass are realized in the gauge-mediated SUSY-breaking model. However, we also find that 
such a condition spans a very small region in the parameter space. Moreover, we also find 
that it is difficult to explain the present baryon asymmetry for such a parameter region, 
unless one fine-tunes the CP- violating A-terms in the potential. 

The paper is organized as follows. In Sec. HIl we give a brief review of the dynamics 
of the Affleck-Dine baryo/leptogenesis and the properties of the subsequently produced Q- 
balls, particularly paying attention to the decay of Q-balls. In Sec. Illlt we evaluate how 
many GWs are produced at the formation of Q-balls and are diluted during the cosmic 
history, which yield the present amplitude of the GWs from the Q-ball formation. Section 
IIVI is devoted to discussion and conclusions. We concentrate on thermal log type Q-balls 
in the main body of the paper. The cases with other types of Q-balls, in which the zero- 
temperature effect (and the negative thermal log effect) dominate the effective potential are 
discussed in Appendix [Bl 



II. AFFLECK-DINE MECHANISM AND Q-BALLS 

In this section, we briefly review the AD mechanism, and the formation and the fate of the 
nontopological solitons, Q-balls. In supersymmetric theories, there are many flat directions 
along which the scalar potentials become flat in the global SUSY limit, and hence such 
scalar fields can easily acquire large field values. The flat directions in the MSSM consist 
of squarks, sleptons, and Higgs fields, and are parameterized by composite gauge-invariant 
monomial operators such as udd and LLe Thus, flat directions can carry baryon {B) 



or lepton (L) charges in general. This is the reason why the AD mechanism can be one of 
the powerful models of baryogenesis. We will see its details in the following subsection. 

The dynamics of a flat direction can be expressed in terms of a scalar field $ (the AD field). 
We consider the dynamics of a flat direction $, taking into account the finite temperature 
effects since these effects play important roles in the formation and the evolution of Q-balls. 
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A. Affleck-Dine mechanism 



The scalar potential vanishes along flat directions in the global SUSY limit. However, 
it is lifted by the SUSY-breaking effects, which depend on the SUSY-breaking mechanism. 
Moreover it is also lifted by a nonrenormalizable operator in the superpotential. 



W 



(1) 



where M is a cutoff scale for an interaction and n is an integer, which depends on a fiat 
direction $. For example, n = 6 if $ parameterizes the udd fiat direction since {uddY is the 
lowest order nonrenormalizable operator. 

In the gravity or anomaly mediated SUSY-breaking scenario, the fiat directions acquire 



their masses 



i 
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where m<^ ~ TeV is the soft SUSY-breaking mass, Mq = I/vSttG is the reduced Planck 
mass, and K is a. numerical coefficient coming from one-loop corrections. The sign of K 
depends on the details of the loop effects: K can be positive if the top quark loop effects are 
the dominant contribution, which is realized when the top Yukawa coupling is order unity. 
On the other hand, K is negative with K = (—0.01 ~ —0.1) when the gaugino loop effects 
are dominant 
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As shown later, Q-balls are formed when K is negative. 
In the gauge-mediated SUSY-breaking model, on the other hand, the potential along the 
fiat direction is lifted as 2^, 



2\ 2 



M| 



m »Ms) 



(3) 



Mp/m ^ is the messenger mass. The allowed 
2^. The upper bound comes from 



where Mp is the SUSY-breaking scale and Ms- 
parameter range of Mp is lO'^GeV < Mp < 10^°GeV 
the condition that the gravity effects should not be so strong and the lower bound comes 
from the condition that the SUSY breaking scale should be larger than the electroweak 
scale. The shape of the potential can be understood by noting that for large |$| > M5, the 
supersymmetry-breaking mass terms are suppressed by a factor of M|/|$p and hence the 
scalar potential becomes constant at |$| > Ms- In addition to Vgaugc? there also exists the 
gravity effects. 



grav2 



3/2 



1 + i^log 



|<|)|2 

Ml 



l$l 



(4) 



where 7723/2 is the gravitino mass ranging from 1 eV to 10 GeV. This term typically dominates 
over V^auge at large field values. 

The scalar potential also contains the contribution from nonrenormalizable operators, 
called y4-term Va and F-term Vp: 



^A = ^^*" + H.c., 



|^|2n-2 
M2"-6 ' 



(5) 
(6) 
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where complex parameter and its absolute magnitude is less than order unity. A- 

term arises from the gravitational interaction between the nonrenormalizable operator and 
the SUSY-breaking sector. 

During both inflation and the inflaton oscillation dominated era, the AD field receives 
the Hubble induced mass term coming from the gravitational interaction between the AD 
field and the inflaton : 

VuM = -CHH^m\ (7) 

where ch is a positive constant of order unity and H is the Hubble parameter. The balance 
between F-term and the (negative) Hubble induced mass term determines the initial value 
of the AD field, but these terms are irrelevant for the subsequent dynamics. 

In addition to the SUSY-breaking effects and the nonrenormalizable operators, there are 
other contributions to the scalar potential. When the thermal plasma exists, the AD field 
receives the finite temperature effects given by j23 |. 



h^T^l^l^ (/i|$|<r), 

|<|)|2 

where h is the Yukawa or the gauge coupling constant for the corresponding AD field, T 
is the temperature of the thermal plasma, c is a numerical constant of order unity, and 
ttg = g'^/4:'7i represents the gauge coupling constant. The sign of c depends on the AD field 
and we assume it to be positive henceforth. The upper term in the right-hand side of Eq. 
dH]) represents the thermal mass from the thermal plasma and the lower one represents the 
two-loop finite temperature effects coming from the running of the gauge coupling g{T) with 
the nonzero field value of the AD field. Note that the thermal plasma exists even before 
the reheating from the inflaton decay. This is because the partial decay of the inflaton 
before its complete decay generates the thermal plasma as a subdominant component of the 
Universe. During the inflaton oscillation dominated era, the temperature of the Universe 
can be expressed as j23| 

T^Af (^^y''(oca-/«), (9) 

where A = 7r'^g^/90 with being the effective relativistic degrees of freedom. The subscript 
"T" indicates that the parameter is evaluated at the scale higher than 1 TeV. Here we have 
assumed that reheating from inflaton decay takes place at the temperature higher than 1 
TeV. Note that, in the context of MSSM, at the energy scale above ~ 1 TeV, (7* ^ 220, at 
the energy scale 100 MeV ~ 1 TeV, g^ ~ 100, at the energy scale 0.1 MeV ~ 100 MeV, 
g^ ~ 10 and at the energy scale below ~ 0.1 MeV, g^: ~ 4. The temperatures at the onset of 
the AD-field oscillation, Q-ball formation and domination, and reheating, which are given 
in Sec. Ill C\ are higher than the electroweak scale when the amplitude of the GWs is large 
enough. Thus, we neglect the time variation of the relativistic degrees of freedom before 
Q-ball formation and reheating and we express A at those epochs as At henceforth. On 
the other hand, the temperatures at Q-ball domination and Q-ball decay can be less than 1 
TeV and hence we express A at that time as A^om and A^ec, respectively. 

Now we consider the dynamics of the AD field. The total effective potential of the AD 
field V is then given by ^ = l^rav/gaugc + (Vgrav2) + Va + Vp + Khcrmai + ^HM, and the AD 
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field obeys the equation of motion, 

dV 

$ + 3/7$ + — - = 0, (10) 

wliere tlie dot denotes tlie derivative witli respect to tlie cosmic time, t. Tlie AD field quickly 
settles down to the potential minimum |<l>| ^ (iirM"'~^)^/*^"~^\ which is determined by the 
balance between Vyoa and Vp. 

VwiA decreases after inflation in response to the decrease of the Hubble parameter and 
disappears after the reheating of the Universe. Then, the AD field (more precisely its radial 
component) begins to oscillate around the origin when 

^osc = y'\^). (11) 

where the dash denotes the derivative with respect to = a/2|$|. Hereafter the subscript 

"osc" indicates that the parameter or the variable is evaluated at the beginning of the 

oscillation of the AD field. The field value at which the AD field begins to oscillate is given 
by0osc^(i^oscM"-3)V(-2). 

The potential of the AD field also contains a phase dependent term, that is, A-term, which 
rotates the AD field unless the phase component of the AD field, Q = arg[$], accidentally sits 
on the potential valley of the A-teim Eq. ([5]). As a consequence, the orbit of the AD field in 
the phase space becomes elliptical and its ellipticity is estimated to be e ^ (imnT'3/2 /V" {(posc) ■ 

If the AD field carries baryon or lepton charge /3c, the angular momentum of the motion 
in the complex plane of the AD field represents the baryon or lepton number density given 
by 

nij(tosc) = i(3ci^*^ - ^*^) - (3camm3/2(j)l,^, (12) 
which implies that baryon or lepton asymmetry is generated in the Universe. 



B. Q-ball formation 

Next we consider the Q-ball formation associated with the AD mechanism. Fluctuations 
around the homogeneous mode feel spatial instabilities and grow nonlinearly during the 
oscillation of the AD field and eventually form clumpy objects, Q-balls, if V{(p)/(j)'^ has 
a global minimum at = ^min 7^ [l^. This is because the pressure of the AD field is 
negative for such a potential [l9|. From Eqs. ([2]), ([3]), and ([H]), the condition can be realized 
when V'gi.av(2) with negative K, V^auge(0 > Ms), or thermal (0 > hT) dominates the potential 
energy. In fact, it is confirmed numerically that fluctuations develop and go nonlinear to 



form Q-balls [20|, |2l|, |23|. Their stabilities are guaranteed by global U{1) charge, that is, 
baryon or lepton charge in our case. In this subsection we briefly review the amplification 
of the fluctuations of the AD field and the properties of the subsequently produced Q- 
balls. Here we concentrate on the case when the dynamics of the AD field is driven by the 
thermal logarithmic potential Eq. ([H]). In Appendix [Bl we will comment on the cases when 
the effective potential is dominated by zero-temperature terms (and a negative thermal log 
term) . 

First, let us examine the growing of the fluctuations of the AD field using the linear 
perturbation analysis. We write the AD field $ as 

$ = ^^^^ (13) 
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and decompose the radial (0) and the phase (6) components into their homogeneous parts 
and perturbations, 



{x,t) = e{t) + 6e{x,t). 



(14) 
(15) 



Once the baryon or the lepton number is fixed, the phase dependent term in the potential 
is irrelevant. Then, neglecting such terms in the potential, the equations of motion in the 
fiat FRW universe read llSlI 



3H(f) 



+ V' = 0, 



e + 3He + 2-e = o, 

{26(1)66 + 6'^6(j)) + V"6(j) = 0, 



66 + 3H66 + 2 I ^50 + ^66 - ^6(j) ] - \a66 = 0. 



(16) 
(17) 

(18) 
(19) 



To find the instability bands, let us write the perturbations as 

50(a;,t) = 50oe^W+*^-*, 
5^(a;,t) = 50oe^«+*^-^, 



S^^\3H 



S 



+ 3HS 



+ Au'S S 



(20) 
(21) 



where 50o and 56^0 are constants. Inserting these forms into Eqs. (fT6|l . (fT7|l . ( ITSjl . and (fT9|) . 
we have the following dispersion relation. 



0, (22) 



where k"^ = and u = 6. The fiuctuations with the momentum k grow exponentially 
when the condition 

— + V"-Lo^< 0, (23) 

is satisfied. Here we have assumed the infiaton oscillation dominated era. The wavenumber 
of the maximal growth mode is given by 



^2 



16a;' 



;{7u''-QV"u'-V 



r/l2\ 



and the fastest growing rate /3gr is given by 
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co^ - V" 



< 



(24) 



(25) 



Finally, such growing fiuctuations go nonlinear to form Q-balls. The Hubble parameter 
at the Q-ball formation may be expressed as 



1 



a 



(26) 



where a > 1 is a numerical factor that represents the duration of the Q-ball formation, and 



a ~ (9(10) for the case of a mass term with negative K [18| and a ~ OiX) for the case of a 



logarithmic term [23 



Next, we investigate the properties of Q-balls. The properties of Q-balls are different for 
different types, and here we concentrate on the case when the oscillation of the AD field is 
driven by the thermal logarithmic potential, 

^thermal ^ «jTMog (^^^ . (27) 

For the gravity or anomaly mediated SUSY-breaking mechanism, this situation can realized 
when 

ajTl » mj0L, (28) 

and for the gauge-mediated SUSY-breaking mechanism, the thermal logarithmic potential 
dominates when 

0osc} for 0OSC > Ms, 

Vosc»|^2^2^^ for 0osc<M5. ^^^^ 

The angular velocity of the homogeneous Affleck-Dine field, cj, is determined by the mass 
scale around (/)osc and is estimated to be \/2a gT"^^^ / (p^sc, which yields from Eq. f l2^ and Eq. 

m 

a' -2 0L ^"-v^0osc- ^ ^ 

The Hubble parameter at the Q-ball formation, Eq. f l26|) . becomes 

1 a 

~ L'fi^^ (31) 

a 0OSC 

and from Eq. ([9]) and Eq. ffTTj) . the temperature at the Q-ball formation is 

T, ~ a-i/^Tosc. (32) 

Almost all the baryon/lepton charges produced by the AD mechanism are absorbed into 
Q-balls. Then, the baryon/lepton charge stored in a Q-ball is estimated as 

e-'-^fn^) . (33) 

\ag ioSC / 

where e is related to the ellipticity of the orbit of the AD field, e, as e = e for e > 0.06 and 
e = 0.06 for e < 0.06, and /3 ^ 6 x 10~^ j23|. The numerical factor /3 represents the dilution 
due to the cosmic expansion during the Q-ball formation. Here one should notice that the 
total charge trapped by a Q-ball is proportional to the ellipticity e for e > 0.06. On the 
other hand, it saturates for e < 0.06, since both negatively charged Q-balls and positively 
charged ones are formed. Other properties of Q-balls are given by. 



47rV2 



^ ay^ToscQ^/^ Eq ^ ^al/'ToscQ'/\ (34) 
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where R is the radius of Q-balls, 0q is the value of at the center of the Q-balls, and Eq is 
the energy per a Q-ball. From Eq. (134|) . the average energy density of Q-balls is estimated 
as 

Pq ^ E^nc ^ ' (35) 

where nq ^ (/cmax/o)^ is the number density of Q-balls and 77 is a numerical factor of order 
unity. From Eqs. (!26l) and (135!) . the corresponding density parameter is given by 

''«-3M^if2- 3^ ^mJ • ^^^^ 



C. The fate of Q-balls 



In this subsection, we examine the fate of Q-balls, which is important for estimating the 
present amplitude of the GWs from the Q-ball formation. Since the dominant contribu- 
tion to the potential depends on the temperature which changes with the cosmic time, the 
properties of Q-balls change as well. Moreover, for some temperatures, other contributions 
can dominate the thermal logarithmic contribution in the potential, which implies that the 
properties of Q-balls may drastically change and Q-balls may disappear if the dominant 
contribution of the potential does not allow a Q-ball solution. Finally Q-balls can decay 
from their surfaces if the energy per charge is smaller than the nucleon mass or the neutrino 
mass. In the following, we investigate the fate of Q-balls in detail. 



1. Transformation of Q -balls 

As mentioned above, the potential for the AD field has temperature dependence and so 
do the properties of Q-balls. The AD-field value at the center of the Q-balls and the energy 
per a Q-ball are given by 

0g(T) ^ «J/2TgV4 ^ (,-^)l/4^_^ ^ (37) 



- osc 



(38) 



Here we have used Eq. fl33|) . We can see that 4>q(T) and Eq(T) are proportional to the 
temperature of the thermal plasma and hence decreases with the cosmic time. Therefore, 
the average energy density of Q-balls is proportional to Ta~^. Since T oc a~^^^ in the 
inflaton oscillation dominated era, so the energy density of Q-balls decreases in proportion 
to oc T^: 

/ osc 

After reheating, it decreases in proportion to Ta^^ oc T^. Thus, as long as Vthcrmai dominates 
the potential of Q-balls, Q-balls never dominate the energy density of the Universe. 
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a. Gravity or anomaly mediated SUSY-hreaking model 



While Vthermai IS propoitional to T^, Vg, 



-m'^(f)Q(T) is proportional to T^. 



grav — 2 

Thus, in the case of the gravity or anomaly mediated SUSY-breaking model, using Eq. ([3 
V^rav dominates thermal at the temperature below Tc given by 



Tr 



1/4 



T 



(40) 



Once V^rav becomes dominant in the potential of the AD field, the properties of Q-balls 
are changed. Figure [T] shows the schematic form of this transition. As the temperature 
decreases, Vthermai and 0q also decrease. At the temperature Tc, the dominant contribution 
to the potential is changed. In the gravity or anomaly mediated SUSY-breaking model with 



Vgrav (0) 








FIG. 1: The scalar potential V^rav + ^thermal is shown. The horizontal axis is the amplitude of the 
AD field (j) = v^l'^l- The filled circle represents the position of the AD field. Initially its rotation 
is driven by Vthermai at T = rosc(> Tc). At the temperature T = Tc, Vgrav dominates over Vthermai > 
and the properties of Q-balls are changed. 

negative K, the potential still allows Q-balls but their properties are changed as follows: 

Eq ^ \m^Q. (41) 

Here, we have used the fact that Q is conserved through the transition. One should notice 
that such Q-balls behave like matter because Eq does not decrease after the transition. In 
the gravity or anomaly mediated SUSY-breaking model with positive on the other hand, 
the Q-ball solution no longer exists, which makes Q-balls unstable. Therefore the AD field 
would be almost homogeneous. The energy density of such an AD field also behaves like 
matter. However, there is an important difference between the cases with positive K and 




10 



negative K. While Q-balls can decay only from their surfaces, the almost homogeneous 
AD field can decay over the whole space. Thus, Q-balls can survive longer and dilute the 
produced GWs further. 

h. Gauge-mediated SUSY-breaking model 

In the gauge-mediated SUSY-breaking model, the situations are rather complicated. 
While thermal is proportioual to and V^rav2 is proportional to T^, V^auge is independent 
of T. Thus, there are three possibilities depending on which term dominates the potential 
next. 

• Case A: Vg^ngc driven Q-ball transformation. 

One is that Vgauge dominates Vthcrmai at the critical temperature given by 

= a-^/^Mp, (42) 

so that the type of the Q-ball is changed into the gauge mediation type. The properties 
of this type of Q-balls are given by 

R ~ —= , u ' 



^ M^gl/^ Eq ^ ^MpQ^'\ (43) 

Note the fact that Q is conserved through the transition. In this case Vgrav2 never 
becomes the dominant contribution in the potential of the AD field. 

Case B : V"grav(-f^ < 0) driven Q-ball transformation. 

Another is that V^rav2 with negative K dominates at the critical temperature given by 



Tc = ^3/2 — , (44) 



SO that the type of the Q-ball is changed into the new type. The properties of this type 
of Q-balls are given by 

2 2 

^ - \KV ^-^3/2, 

''''3/2 I-"- I 

0Q ^ i^—ij m3/2Q'/^ Eq ^ -my^Q. (45) 

Here, Q is also conserved at the transition. In this case, V^auge never becomes the 
dominant contribution in the potential of the AD field. 

Case C : V^rav(-^ > 0) driven Q-ball transformation. 

The other one is that V'grav2 with positive K dominates first at the critical temperature 
given in Eq. (jH]). This potential does not allow a Q-ball solution and hence 
the almost homogeneous AD field is recovered. After this transition, the AD field 
decreases as oc because of the cosmic expansion. Thus, V^auge dominates Vgrav2 at 
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= 0eq = Mp/'m^i2i which implies that Q-balls are formed again. The properties of 
this type of Q-balls are those of the delayed type Q-balls and are given by 



R 



Q 



1/4 



y/2MF 



(v^m3/2)"\ 



V2: 



7rm3/2, 



Q 



(46) 



In this case, Q is given by 



Q 



■'eq 



'If 



Mf 



(47) 



In all cases, Q-balls behave like matter because Eq does not decrease after the transition 
and hence Q-balls may dominate the energy density of the Universe. Figure [2] shows the 
schematic form of the transition of each case. As the temperature decreases, thermal and (pq 
also decrease. The dominant contribution to the potential is changed at the temperature 
Tc^^\ and hence the type of Q-balls are changed as well. 



14™. (0) 




Case A 



^ 



V,rs. (0) 




Case B 



^0 



(0) 




^0 



FIG. 2: The scalar potential in the gauge-mediated SUSY-breaking model in each case is shown. 
The horizontal axis is the amplitude of the AD field. The filled circle represents the position of 
the AD field. Initially its rotation is driven by I4hermai at T = Tqsc and dominant potential terms 
change at T = Tc^^^- There are three possibilities depending on which term dominates afterwards. 
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2. Decay of Q-halls and AD field 



Q-balls can decay into light fermions if the decay processes are kinematically allowed. 
However, in their interiors the Pauli exclusion principle forbids their decays into fermions 
(28| . Therefore Q-balls can decay only from their surfaces. This sets the upper bound on 
the decay rate of Q-balls, 



dQ 



dt 



< 



487r ' 



(48) 



In fact, it is almost saturated for the cases we are interested in 
the decay rate of Q-balls as 

^ _ 1 dQ _ u^R^ 



28(1 . Then, we can express 



dec 



(49) 



487rQ' 

and the Hubble parameter at the Q-ball decay is given by i^dec = Tdec^- The decay rate 
of the Q-balls is different according to their types. In the case of the gravity or anomaly 
mediated SUSY-breaking models with negative from Eq. (14T1) . the decay rate of Q-balls 
can be expressed as 



r =1^ 

J- dec — ^ 



1 



Q dt Q24n\K\' 

As mentioned above, in the gauge-mediated SUSY-breaking case, there are three possibilities 



for the fate of Q-balls. From Eqs. (Il3 



6]), and fl48|l . the decay rate can be expressed 



as 



dec 



1 dQ 

qIE 



1 Tl'^Mp 

Q2Atx\K\ 
, 24\/2M> 



for Case A, 
for Case B, 

for Case C. 



(51) 



The decay temperature should be higher than 1 MeV for the successful big bang nucleosyn- 
thesis (BBN), which constrains the parameters of Q-balls, such as (pose and Tr, if Q-balls 
dominate the energy density of the Universe.^ The constraint in each case will be given in 
the next section. 

Note that Q-balls can evaporate when there are thermal plasmas, as discussed in Ref. 
(29| . Though Q-balls can evaporate away before their decays, this takes place only when 
Q is small enough, which is unfavorable for our purpose because the initial amplitudes of 
the produced GWs would also be small. In Appendix A, we will give more quantitative 
discussion on the charge evaporation from Q-balls. 



^ If there are light bosons and the AD field can decay into them by two-loop interactions, the decay width 
would be enhanced by a factor < 10"^. However, such bosons do not exist in the context of MSSM 
because scalar fields other than the flat direction acquire large masses in general, and hence we do not 
consider such decay processes. 

^ More precisely, if Q-balls contribute to more than about 1% of the energy density of the Universe, their 
decay temperature should be higher than 1 MeV for the successful BBN. But, this constraint is not so 
important for our estimate, and we ignore such small difference. 
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On the other hand, in the case of the gravity or anomaly mediated SUSY breaking models 
with positive Q-balls vanish at T = Tc and hence the almost homogeneous AD field is 
recovered, whose amplitude decreases as oc if due to the cosmic expansion. As long as 
h(f)Q > m^, the fields coupled with the AD field acquire large masses so that the AD field 
cannot decay into them. When h(j)Q ^ m^, the decay into the light fermions is allowed 
kinematically. The decay rate is given by 

Tdec ^ -^m^f- (52) 

In the case where the AD field dominates the Universe at that time, the Hubble parameter 
is given by 

- Jk- '''' 

Therefore, if /i > {tti^/MgY/^ ~ 10"^, then, Tdec > -f^dec, and the AD field decays into light 
fermions quickly at that time. The decay temperature T^^c is estimated as 

^ - >m^. (54) 



Thus, the AD field can decay before the BBN and the electroweak symmetry breaking 
Before closing this section, we comment on the decay rate of Q-balls adopted in Ref. [1 



The authors of Ref. [IJ] considered two possibilities of the Q-ball decay to enhance its decay 
rate. One is to introduce higher-dimensional operators, which preserve the SUSY but violate 
the baryon numbers. However, if such operators provide the decays into fermions, the decay 
rate is strongly constrained by the Pauli blocking as mentioned before. We also would like 
to point out that it is difficult for such operators to provide the decays into bosons in the 
context of MSSM, because such bosons acquires large masses through the Yukawa couplings 
so that the decays into them are kinematically prohibited. Furthermore, even if we find 
such bosons, the decay rate estimated in Ref. cannot be applied to Q-balls, as pointed 
out by the authors themselves. The other is the semiclassical decay through the A-terms, 



which was investigated in detail in Ref. [30|. In Ref. [30[, it is found that the instability 
due to the A-term causes the Q-ball decay in the Minkowski background. Then the authors 
of Ref. [11 simply apply the expression of the decay rate to the case of the expanding 



Universe, which is not justified, and conclude that Q-balls quickly decay via this process. In 
the expanding Universe, however, this process is not effective because the amplitude of the 
AD field decreases and the A-teim. becomes soon negligible due to the cosmic expansion. 



In fact, the authors of Ref. [30| conclude that Q-balls cannot decay via this process in the 



expanding Universe. Therefore, in this paper, we use the conservative estimate of the decay 



rate of Q-balls, different from Ref. [14 



III. GRAVITATIONAL WAVES FROM Q-BALLS 

In this section we investigate the features of the GWs produced from the Q-ball formation 
and discuss the prospect for the detection of such GWs. First we estimate the amplitudes 
and the typical frequencies of the GWs at the Q-ball formation. Then, we evaluate the 
dilution of the GWs during the subsequent cosmic expansion. As we have seen in the 
Sec. Ill CI Q-balls driven by the thermal logarithmic potential never dominate the energy 
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density of the Universe until their transformation. For this reason, we concentrate on the 
situation where Q-balls are formed when the thermal logarithmic potential is the dominant 
contribution to the potential of the AD field. The other cases where Q-balls are produced 
for the zero-temperature potential are discussed in Appendix B. 



A. Generation of GWs 

In this subsection, we study the generation of the GWs associated with the fragmentation 
of the AD field and estimate the amplitudes and frequencies of the GWs. The process of 
the fragmentation of the AD field is inhomogeneous and nonspherical so that GWs may be 
emitted at the Q-ball formation. 

We evaluate the initial amplitudes and frequencies of the GWs from the Q-ball formation, 
by using the equations for the transverse-traceless (TT) component of the metric perturba- 
tions, following Refs. lH, l3| . By perturbing the Einstein equation, we obtain the equation 



for the TT component of the metric perturbations, 

hj^{x, t) + 3HhJ^{x, t) - -^hj^ix, t) = lenGUj^ix, t), (55) 

where hj^{x,t) is the TT component of the metric perturbation and I{J^{x,t) is the TT 
component of the energy-momentum tensor of the AD field. Instead of using the above Eq. 
(!55|) . it is easier to use the following equations: 

Uij{k, t) + SHiiijik, t) + —Uij{k, t) = lQnGTij{k, t). (56) 

Qi 

Here Uij{k,t) and Tij{k,t) satisfies the relations 

^ij i^i^) ~ ^ij,mnik)Umn{k,t), (^'^) 
nJ.T(fc, t) = A,j^mnik)Tmnik, t) , (58) 

T,j{x,t) = ^di<l){x,t)dj(P{x,t), (59) 
where the projection tensor Ajj^^n is defined by 

^ij,mn{k) = ^Pim{k)Pjnik) — — Pjj (fc) (fe)^ , (^O) 

Pij{k) = 5ij - kikj, (61) 

with ki = ki/\k\. hj^{k,t) and n^'''(fc,t) are Fourier transforms of hjj^{x,t) and Ujj^{x,t), 
respectively. Since Eq. ( I56l) contains unphysical (gauge) degrees of freedom, we need to 
follow the time evolution of the variable /i^-^ instead of Uij, strictly speaking. However, in 
the absence of spherical symmetry and homogeneity, it is sufficient to approximate hj^ by 



Uij 
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The energy density of the GWs is given by 

1 

327rGL3 



PGW = TjTT^TT^ / d^xhj;^ix, t)h^^{x, t), (62) 
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which can be approximated as 



where V = is the volume of the space, and we have used the fact that the process of 
fragmentation is nonspherical and \hj^\ ^ \uij\ in this situation. 

We now estimate the energy density of the GWs produced by the fragmentation of the 
AD condensate. For this purpose, we first evaluate the energy-momentum tensor of the AD 
field Tij. The maximal growth mode of the fluctuation of the AD condensate evolves as 

5(t){x, t) = 50oe^^'*+^^"^^''-^, (64) 

where 6(j)Q is the initial value of the field perturbation. Then, Tij can be estimated as 
Tij ~ —k'^^^6(f)'^/3a^ after angle average, and Eq. fISBl) becomes 

4^.2 2 ^2 

W.i(2fcmax, t) + -^Uij{2k^,^, t) ^ -^^502^2/3,.*^ (g5) 

G 

where the Hubble friction term is neglected since k^sa^/a is larger than H. Since /3gr is 
smaller than Ajmax/o (Eq- (I25D), "we may use the second term in the left-hand side in Eq. 
( l65l) to estimate Uij to give 

After a time interval At ^ ln(0Q/50o)//3gr, the fluctuations of the fiat direction reach 5(j) ~ 
0Q, which implies that Q-balls are formed. Therefore, at the last stage of the formation of 
Q-balls, pgw reaches the maximal value, which is estimated using Eqs. (!63|) . (164|) and ( 166|1 

as 



M2. . Pi 



Pgw ^ —Ui.iiij ~ 77^0Q. (67) 



The typical frequency of the GWs from the Q-ball formation is given by the wave number of 
the maximal growing mode, f\ ~ A;max/(7ra). The density parameter of the GWs, ^Gwif*)-, 
is then given by 

"cw(/.) - (68) 

where if* is the Hubble parameter at the Q-ball formation. The present density parameter 
and the frequency of the GWs from Q-balls are then given by 

(^l)' (I) \ (69) 

^» - f' (S^) ■ f^") 

Therefore, in addition to the redshift factor, the present density parameter contains the 
dilution factor due to the matter/Q-balls domination effects. 
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In the case when the oscillation of the AD field is driven by the thermal logarithmic 
potential, from Eqs. ( 130|) . (13T]) . (133|) and ( 134|) . the density parameter and the frequency of 
the GWs at the Q-ball formation are given by 

- (^^) 

It may be instructive to give alternative derivation of the energy density of the GWs 
using the quadrupole approximation. By taking the moment of inertia as / ~ EqE? ~ 
0QCi;^(/cniax/a)~^ and approximating the time derivative by /3gr ~ k^ax/d ~ u, the quadrupole 
approximation gives the luminosity Lqw — / Mq and thus energy density liberated in 
GWs during the Q-ball formation (the interval At ~ (3~^) as 

where we have used (130|1 and (l34l) and nq ^ (A^max/o)^ ~ Apart from the numerical 
factor, Eq. (173|) nicely coincides with Eq. (|6j 



B. Cosmic history 

In order to evaluate the present amount and typical frequency of the GWs from the Q-ball 
formation, we need to take into account of the cosmic history after the Q-ball formation. 
In the case where Q-balls are formed when the thermal logarithmic term is the dominant 
contribution to the potential of the AD field, there are four possibilities of the cosmic history. 
They are classified according to the following two conditions. The first condition is whether 
the Q-ball dominated era exists or not, and it is characterized by i^dom and -f^dec- -f^dom is the 
Hubble parameter when the Q-balls (or the AD field) would dominate the energy density of 
the Universe if such an epoch would exist, i^dec is the Hubble parameter when the Q-ball (or 
the AD field) decays. The second condition is whether the reheating after the infiaton decay 
occurs before the transformation of the Q-balls due to the change of the effective potential 
or it occurs after. This condition is characterized by Tr and T^. Here Tr is the reheating 
temperature of the inflaton and Tc is temperature at which the zero-temperature potential 
dominates over the thermal logarithmic potential and hence the properties of the Q-balls 
are changed. The four cases are shown in Table [B The time evolution of the energy density 
of each component of the Universe in each case is shown in Fig. [31 The dilution factor of 
the GWs from Q-ball formation is different in each case. 



Case 1 


Reheating =^ Q-ball transformation =^ Q-ball/AD-field domination =^ Q-ball/AD-field decay 


Case 2 


Reheating =^ Q-ball transformation=^ Q-ball/AD-field decay (No Q-ball/AD-field domination) 


Case 3 


Q-ball transformation =^> Reheating =^> Q-ball/AD-field domination =^ Q-ball/AD-field decay 


Case 4 


Q-ball transformation =^ Reheating =^ Q-ball/AD-field decay (No Q-ball/AD-field domination) 


TABLE I: The possible cosmic histories. 
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9h 3c 



a 



Case 1 




3c 3de 

Case 2 



3 




3(; ^dom 5de 

Case 3 




3c 3r 3ie 

Case 4 



FIG. 3: The schematic time evolution of the energy density of each component (inflaton, radiations 
from the inflaton decay, Q-balls, radiations from the Q-ball decay, GWs) for each case in Table lU is 
shown. The vertical axis represents the energy density and the horizontal axis represents the scale 
factor a. The dilution factor of the GWs from the Q-ball formation is different in each case. Even 
during the inflaton oscillation dominated era, the energy density of the thermal log type Q-balls 
changes as the temperature of the thermal plasma generated by partial inflaton decay decreases. 



We note that the following relations for if* and Tqsc, which are derived using Eqs. 
f lTT]) . (P7|) and ([SI]), are useful in the subsequent discussion, 



aA. 



1/2 



/,2 
osc 



and 



1/4 



agMc 



1/2 



R- 



(74) 



C. Gravity or anomaly mediated SUSY-breaking model 



In this subsection, we evaluate the present amount and the typical frequency of the 
GWs from the Q-ball formation in the case where the Q-balls are formed when the thermal 
logarithmic term is the dominant contribution to the potential of the AD field in the gravity 
or anomaly mediated SUSY-breaking model. In this case, from Eqs. (128!) and (1741) . we have 
the lower bound on the reheating temperature 



^1/4 



an 



M, 



G 



1/2 



^osc — '^R- 



(75) 
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This bound strongly constrains the frequencies of the produced GWs. 



1. Gravity or anomaly mediated SUSY-breaking model with positive K 

First, we consider the gravity or anomaly mediated SUSY-breaking model with positive 
K. In this case, when the zero-temperature potential dominates the effective potential, the 
Q-balls become unstable and hence decay. Therefore the almost homogeneous AD field is 
recovered at the critical temperature Tc, Eq. PUI) . as discussed in Sec. Ill C. 

First of all, we investigate the condition Tr > T^. in detail, which discriminates Cases 1 
and 2 (Tr > T^) from Cases 3 and 4 (Tr < Tc). Inserting the expression of T^sc Eq. 
into Eq. (jlD]) yields 

Tc = (76) 

which shows that T/j > is equivalent to 

.1/8 m^^^S^^^ 
ag' MJ 

Next, in order to examine the condition of the existence of the Q-ball dominated era, we 
evaluate the Hubble parameter at the would-be Q-ball domination, i^dom- In fact, Hdom in 
Case 1 and that of Case 3 coincide with each other. This can be understood as follows. 
Qq oc T during the inflaton oscillation dominated era with thermal log type Q-balls and 
Qq oc during the radiation dominated era with nonthermal type Q-balls, while Qq 
is constant both during the inflaton oscillation dominated era with zero-temperature type 
Q-balls and during the radiation dominated era with thermal log type Q-balls. Thus, both 
for Case 1 and Case 3, the following relation is satisfied: 



= = 1. (78) 

dom 

Therefore we have 



rp2rp2 

TT ^dom iJ-' c o*2 (>7n\ 

and from Eqs. ([32]), (jSS]), dZH) and ([76]), we obtain 

= (80) 

Here the subscript "dom" indicates that the parameter or variable is evaluated at the would- 
be Q-ball domination. On the other hand, i^dec is given by Eq. (I53l) . Comparing these 
equations, the condition of the existence of the Q-ball dominated era is equivalent to 

6V^a^/^/zVM;l4/^ _ 

^'^^Mj = ^^- (^^^ 

From Tc in Eq. ( 176]) . i/dom in Eq. ( [80]) . and if dec in Eq. ([53]) . we thus obtain the following 
conditions on Tr for each cosmic history, which are summarized in Table [Tl| 
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Case 


1 




Case 


2 


Tr > max{T^n0osc),Tf (0osc)} 


Case 


3 


Tr < min{r|.i(0osc),T^^(0osc)} 


Case 


4 





TABLE II: The conditions of four cases of the cosmic history. 



We then calculate the present density parameters of the GWs, f^QW' ^^"^ their frequencies, 
/o, for each case. 



D 



a* a. 



eq 



gr Go 

Hr 



2/3 



(Case 1,3), 
(Case 2,4), 

2/3 



dec 



TJ \ 2/3 
Hr \ a 



dom 



eg 
ao 



(82) 



(Case 1,3), 
(Case 2,4). 



(83) 



Here D represents the dilution factor of the GWs due to the matter /Q-balls domination. 
The redshift at the Q'ball formation is given by 



a* _ 
ao 



At 
Ao 
At 



1/3 



1/3 



4 \ 1/12 ^ 



To 
Tr 



^dcc 

Hr 



Tr \H, 



2/3 



H, 



dec 



' dom 



1/6 



Hr 

H^ 



2/3 



(Case 1, 3), 
(Case 2,4), 



(84) 



where the subscript "0" indicates that the parameter is evaluated at present. Here we have 
neglected the difference between and g^:s, where g^s is the relativistic degrees of freedom 
for the entropy density, aeq and ao are the scale factors at the matter-radiation equality and 
at present, respectively, and so that ao/a^^ ~ 3200. From Eqs. ([53]), (JH]), dZZ]), (EH), flHOjl . 
f l83|) . and fl8^ . we thus obtain the present density parameter of the GWs, I^qw? 



"GW — 



18 • 2V3 a^^A^ilhy^ 



-1/3 



M, 



G 



16/3 



''eq 



ao 



and the present frequency, /o, as 



/o 



R 



Mr 



4/3 



a, 



eg 
ao 



(Case 1,3), 
(Case 2,4), 



2^/1% a 



271^/' 



To 



/,2 
osc 







(Case 1,3), 
(Case 2,4). 



(85) 



(86) 



20 



In Fig. m we show the contour plot of the present density parameter of the GWs, f^cw) 
and their frequency, /q. Although the amplitudes of the GWs can be large, their typical 
frequencies (> 10'^~^ Hz at least) are too high to be detected by the future experiments, since 
the sensitivity ranges of future detectors are ^ 10~^ fo — 10^'"^ Hz for advanced 
LIGO, n^y^ > 10-1^ at fo ~ 10-3~-2 for LISA, and fi^^ ^ 10"^^ at fo ~ IQ-^ ~ 10 
Hz for DECIGO. Large f^GW requires high reheating temperature and high decay rate, but 
that results in higher frequency of the GWs in turn. 




FIG. 4: The contour plot of the present density parameter of the GWs, J^qyv ~ 10~^^, 10"^'* and 
-|^Q-i6 lines) and its frequency, fo = 10^ Hz, 10^ Hz, and 10^ Hz (green lines) is shown. Dashed 
lines represent and T^, which distinguish Cases 1-4 (Table [Tl]). In the region surrounded by 
T = and <^osc = Mq, there is the AD-field dominated era. Dotted line represents the boundary 
of the evaporation of Q-balls, Eq. (j88|) . In the shaded region, thermal log type Q-balls are not 
formed, Tr < T^. The parameters are taken to be r/ ~ l,/i ~ 0.1, a ~ 1, ~ 0.1, and ~ 10^ 
GeV, and the effective relativistic degrees of freedom of MSSM are assumed for At, ^donu ^dcc and 

In the above argument, we do not consider the possibility that the Q-balls evaporate 
out via diffusion processes before their decay. However, that could be possible only when 
charge Q stored in a Q-ball is very small. This can be understood as follows. First one 
should notice that from Eqs. fl33|) . ([71]), and f l76|) . the typical charge of the Q-ball and the 
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transformation temperature are given by 

~ 10-4 '^osc ^ J, _ m^(j)lic , . 

Inserting these values into Eq. (1A5I) in Appendix A yields 

AT-i2/3 Ti^7/15 / rp \ 2/3 —2/15 

< 10«/-Zi^ . lO-GeV (^) (^)- .Cf. (88) 



In Fig. m this constraint is represented by dotted line above which Q-balls evaporate out 
before their transformation. From Fig, HI we find that the evaporation is effective at /o > 10^ 
Hz for i^QYV > 10~^^. Therefore the charge evaporation from the Q-balls does not change 
our conclusion. 



2. Gravity or anomaly mediated SUSY-breaking model with negative K 



Next we discuss the gravity or anomaly mediated SUSY-breaking model with negative K. 
The only difference between the cases with positive and negative K is that Q-balls do not 
disappear even after the critical temperature Tc for negative K. Thus, the cosmic history 
can be classified to four possibilities in the same way as the case with positive K. However, 
the decay time of the Q-balls, -ffdec, is different from that of the AD field for positive K. 
From Eqs. fl33|) . f l50|) . and fl7i|) . it can be estimated as 



H, 



dec 



(89) 



24:7i\K\At 

In contrast to the case with positive K, the decay temperature ( ^dec^ V -^dec ) can be less 
than 1 MeV, which is forbidden by the BBN constraint if Q-balls dominate the Universe 
before their decay. Thus, we have another constraint on the reheating temperature as 

^ 1/4 



Tp > 0.4 X 



dec 



mm 



M, 



lO^GeVy 



Mg = T\ 



c,BBN 
R 



(90) 



for Cases 1 and 3. 

Using Eqs. dH]), ((72]), ([74]), ([80]), (|83]), ([84]) with the decay time Eq. ([89]), the present 
density parameter of the GWs, I^gW' given by 

l4/3,v4 /l#2/3. 



1 

216 
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a'I^A\'LAT\K\^'^ 



,2/3 
^dom^T 
16/3 

•'osc \ 



-Lr "eg 



m 
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,4 ao 

osc 
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eg 

ao 



and the frequency, /o, is given by 

1/6 1/3^2 ^1/4 .1/3 
V "q^dom^O 



fo 



73 m 
2tt \7r) 



a^/^Al" 



.1/12 
^dec 

Tr 



aVl2^Vf^3/4|^|,/6 



-l/2/i?^'^G . 



1/6 , 



,3 

osc 



27rA 



1/6 
T 



.2/3^1/3 



(Case 1, 3), 
(Case 2,4), 

(Case 1, 3), 
(Case 2,4). 



(91) 



(92) 
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From Eqs. (1761) . (IHOi) . and (189|) . the critical temperatures T]^^ and in this case are given 
by 



t: 



cl 
R 



.1/8 m^/^(^^/^ 



3/4 



M, 



1/4 
G 



o lr^l\l/6 3/44I/12.I/3 m^/^^'''/^ 



1/3 4/3 



3/2 



(93) 
(94) 



G 



which characterize all CclSGS cLS summarized in Table [ITl 

Figure E] shows the contour plot of the present density parameter of the GWs, f^GW' ^^"^ 
their frequency, /q. f^cw independent of the reheating temperature, Tr, in Cases 2 and 4 
while it depends on both T/j and 0osc in Cases 1 and 3. /o is sensitive to Tr. We find that it 



GeV 




osc 



GeV 



FIG. 5: The same plot as FigU] but for the case of the gravity-mediated SUSY-breaking model 
with negative K is shown. The shaded region represents the forbidden region where Tr < or 
Tr < T^'^^^. The parameters are taken to be 77 ~ 1, a ~ 1, ~ \K\ ~ 0.1 and ~ lO^GeV, 
and the effective relativistic degrees of freedom of MSSM are assumed for Ay, y4domi ^dcc and Aq. 



is almost impossible to detect the GWs from Q-ball formation even by the future detectors 
in this case, as with the case with positive K. The relatively small decay rate of the Q-balls 
leads to a lower frequency of the GWs, but it also results in the Q-ball domination, which 
dilutes GWs. 
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We comment on the possibility of the Q-ball evaporation. In this case we need to consider 
the Q-ball evaporation both in the thermal log case [Eq. ( lASP ] and in the gravity-mediation 
case [Eq. ( ]A9P ]. The condition for the evaporation of the Q-balls before their decay in 
the case when the thermal logarithmic term dominates the potential is the same as that 
for positive K. Using Eq. (18 7p . the condition for the Q-ball evaporation after the Q-ball 
transformation is 

'^"^^ < ^ ^4/27 - \ io9GeV ) \ ITevJ " ^^^^ ' ^^^> 



As in the case with positive K, we find from Fig. |5] that the evaporation of Q-balls can 
modify the above estimation of and /o only when /o is too large (/o > 10^ Hz). 

Therefore, even if we take the evaporation effects into account, the conclusion is unchanged. 

Summarizing our results in the gravity or anomaly mediated SUSY-breaking models, we 
find that it is possible to generate a large amount of GWs from the Q-ball formation, but 
the present frequency is relatively high. This is mainly because the initial frequencies of 
such GWs are rather large and cannot be redshifted enough. Thus, we conclude that in this 
mediation mechanism, the GWs from the Q-ball formation cannot be detected by the future 
detectors (designed or planned) even if thermal log terms induce the Q-ball formation. 



D. Gauge-mediated SUSY-breaking model 

In this subsection, we study the present properties of the GWs from the Q-ball formation 
in the case where the Q-balls are formed when the thermal logarithmic term is the dominant 
contribution to the potential of the AD field in the gauge-mediated SUSY-breaking model. 

From the condition that the thermal logarithmic term dominates the potential of the AD 
field, Eqs. (I29p . we have the lower bound on the reheating temperature 

'''' 

where we have used the relation Eq. (!74|) . This bound strongly constrains the frequencies 
of the produced GWs. As mentioned before, there are three possibilities of the Q-ball 
transformation. Cases A, B, and C (Fig. [2]), depending on the dominant term (Vgaugc or 
V^rav2) in the potential and on the sign of K. We estimate the present properties of the 
GWs from the Q-ball formation in each case. 



1. Case A: V^augc driven Q-ball transformation 

First we consider Case A, in which the Q-balls are formed by the thermal logarithmic 
potential and their type is changed into the gauge-mediated type at the critical temperature 

A 

From Eqs. flHT]) and f E7|) . the Hubble parameter at the Q-ball decay is given by 
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We have three constraints on the parameters of this modeL The first constraint comes 
from the BBN. For the successful BBN, the Q-balls must decay before BBN and the decay 
temperature should be larger than IMeV if Q-balls dominate the energy density of the 
Universe. This requirement becomes 



.1/4 .1/10 

Tfi>3x 10-^(6/3)^/^ ^ ^""^ 



a. 



,3/2 

yoSC _ rpC,BBN 

■^3/10^1/5- R ■ 



(98) 



The second one is the condition for Vthcrmai to dominate the potential of the AD field Eq. 
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The last one is the condition that Vgaugc > K-mv^ at 

,1/4 



a. 
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Here we have used Eqs. fl37j) and P2|) . 

From Eqs. ([32]), (jSH]), (iJ]), ([71]), and ([79]), obtain the Hubble parameter at the Q-ball 
domination as 



H. 
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^9/241/2 4 1/2 ,5 /l#2 



The present amount of the GWs from the Q-ball formation is then given by 
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and the frequency, /o, is estimated as 
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Here we have used Eqs. ([71]), ([72]), ([71]), ([83]), ([97]), and ffTOT]) . The critical tempera- 

tures, and T^, using Eqs. ( I42l) . (l97l) . and OlOip . can be estimated as 
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which characterize all four cases as summarized in Table [TTl 



25 



The solid (red) line in Fig. E] represents the maximal amplitudes of the GWs, f^Q^'^'', 
from Q-ball formation at the frequency ranging from 1 Hz to 10^ Hz in Case A. When 
/o ^ 10^ Hz, fiQ^*^"^ oc /g''^, which corresponds to the parameter region with Tr = T^^^^ , 



lO^GeV < Mf < 10^°GeV, and 0osc = Mq. When lO^Hz < /o < 10^ Hz, fijj!;;^" oc /( 



which corresponds to the parameter region with T^' < Tr < T^, 0osc = Mq, and 
Mp = 10^ GeV. Note that in both cases, there is the Q-ball dominated era. When /q > 10^ 
Hz, l^Q^^'' is obtained for 0osc = Mq without the Q-ball dominated era. As is seen in Fig. 
El f^GW^^ becomes larger than 10~^^ but at /o > 10^ Hz, which goes outside the sensitivity 
range of the DECIGO or BBO. Therefore, it is almost impossible to detect the GWs from 
Q-ball formation even by the future detectors in this case. 

We comment on the possibility of the Q-ball evaporation. In this case, Q-balls can 
evaporate out before the decay evaluated above and the estimation of f^^w ^^"^ /o can 
break down. However, this is the case only for /q > 10^ Hz. Thus, the above conclusion is 
unchanged by the evaporation effects. 

2. Case B: lgrav(-f^ < 0) driven Q-hall transformation 

Next we consider Case B. In this case, Q-balls are formed by the thermal logarithmic 
potential and change their type into the gravity-mediated type {K < 0) at the critical 
temperature . This case is similar to the case of the gravity- mediated SUSY-breaking 
model with negative K. Thus, we need only to replace by 7713/2 (< rUfp). While the larger 
amount of the GWs is obtained by replacing by 777,3/2 (< m^) of Case 3 in Eq. fl92l) . which 
is favorable for the detection, the present frequency of the GWs is also enhanced, which is 
unfavorable for the detection. 

There are also three constraints on the reheating temperature. The first one is the BBN 
constraint, which is required when Q-balls dominate the energy density of the Universe, and 
is given by 

T« > 2 X m-'A}^'AZ (^) "* (^) (^) Ma . Tf-r (106) 

This constraint can be derived in the same way as deriving Eq. ( 1901) . The second one is the 
condition that the thermal logarithmic potential should dominate the potential of the AD 
field at the Q-ball formation and is given by [Eq. (l96l)]. 



The last one is the condition for Vgrav to dominate Vgaugc at the Q-ball transformation. 



Here we have used Eqs. (1571) and dS]). Since m^i2/MF ~ Mp/Mc ^ 10~^ if there is no large 
hierarchy between the SUSY-breaking sector and the messenger sector, the first constraint 
contradicts with the last one. However, if there is a hierarchy between them, three conditions 
can be made compatible. 




(107) 




(108) 



26 



The dashed (green) hne in Fig. E] represents the maximal amount of the GWs, ^q^'^'', 
from the Q-ball formation at the frequency ranging from 1 Hz to 10^ Hz in Case B. When 
/o ^ 10^ Hz, ^Q^^"" oc /o^^, which corresponds to the parameter region with T/j = T^'^^^, 
m3/2 < 10 GeV, and 0osc = Mg- For lO^Hz ^ < 10' Hz, l]^^'^^ ex ff^\ which 
corresponds to the parameter region with Tr = T'^^^ , 0osc = Af^, and 7723/2 < 10 GeV. Note 
that in both cases, there is the Q-ball dominated era. For /o > 10^ Hz, ^(^^ is obtained 
for 0OSC = Mg without the Q-ball dominated era. As seen in Fig. [6l ^(^^ is larger than 
10"^^ for /o > 10'^ Hz, which again goes outside the sensitivity range of the DECIGO or 
BBO. Therefore, it is almost impossible to detect the GWs from Q-ball formation even by 
the future detectors in this case too. The conclusion is unchanged by the evaporation effects 
for the same reason discussed in Sec. HHP 1[ In this case, the Q-ball evaporation out takes 
place when /o > 10^~^ Hz. 



3. Case C: V^rav2(-^ > 0) driven Q-ball transformation 

Finally, we consider Case C. In this case, Q-balls are formed by the thermal logarithmic 
potential and are transformed into the almost homogeneous oscillating AD field at the 
critical temperature T^. Afterwards, the potential energy of the AD field decreases and 
^augc dominates the potential finally. At that time, Q-balls are formed again. Note that 
the second Q-balls are of the "delayed" type and hence the GWs from the second Q-ball 
formation cannot be detected, as we will see in Appendix. B. 

From Eq. f l5T]) . the Hubble parameter at the Q-ball decay is given by 

i?dec = (109) 

We have again three constraints on the parameters of this model. One comes from the BBN 
constraint, 

m3/2 > 4 X 10~' A'J^fMfM'f ^ m^fl (110) 

which is required when the Q-balls dominate the energy density of the Universe. Another 
is the condition for Vthcrmai to dominate the potential of the AD field, which coincides with 
Eq. dMD , 

,1/4 f my 2 _^e,th 
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Tn>AT\-^] 0osc = Tr. (Ill) 



The last one is the condition that V^rav > V'gaugc at = 0c, which coincides with Eq. f llOSp . 

/1 1/4 i,3/2 

Tn < ^(,/3)V4!!!|^ ^ r^^r, (112) 
From Eqs. fillip and flll2p . we have the inequality 



2 



0OSC > m-'^- (113) 
"^3/2 

If there is no hierarchy between the SUSY-breaking sector and the messenger sector, M|. ~ 
777,3/2^^, the inequality f lllSp cannot be satisfied. Hence we consider the case when there is 
some hierarchy between them so that Eq. (I113p is satisfied. 
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The Hubble parameter at the Q-ball transformation is given by 



dom 
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which is obtained by replacing in Eq. 
from the Q-ball formation is then given by 
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and the frequency, /o, is estimated as 
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Here we have used Eqs. dZI]), ([72]), ([71D, (jHSD, dH, ffTOOD . and ffTTij) . From Eqs. 
and (11141) . the critical temperatures, and T^, are estimated as 
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which characterize all four summarized in Table ITTi 

The dotted (blue) line in Fig. [6] is the maximal amplitudes of the GWs, ^2^^*^^, from the 
Q-ball formation at the frequency ranging from 1 Hz to 10^ Hz in Case C. For /o < 10'^ Hz, 
^Q^^^ oc /o^'''', which corresponds to the parameter region with Mp ~ 10^ GeV, 777,3/2 < 10 



Mg and Tn > T^'*^ for 



GeV and Tr = T^'**^ for 0osc = M^. For 10^ Hz < /o < 10^ Hz, fi^^^^ ex /o, which 
corresponds to the parameter region with Mp ~ 10^ GeV, 
7773/2 = 10 GeV. Note that in both cases, there is the Q-ball dominated era. For /o > 10^ 
Hz, l^Q^'^'^ is obtained for (^osc = Mg and there is no Q-ball dominated era. As seen in Fig. 
il Q^Q^^"" is larger than lO'^^ for /o > 10 Hz, which is on the edge of the DECIGO or BBO 
sensitivity range. Thus, it is difficult but not impossible to detect such GWs by the next 
generation detectors. The parameters that realizes f^Q-vv — 10~^^ at /q ^ 10 Hz are given by 



lO^GeV, 



7«3/2 



lOGeV, 



Mg and Tr 



rpC,th ^ 



10^°GeV. 



(119) 



The conclusion above is unchanged by the evaporation effects for the same reason discussed 
in Sec. IIII D 1[ In this case, the Q-ball evaporation takes place when /q > 10 Hz. One 
should notice that though gravitino does not overdose the Universe because of large entropy 
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10"^ 




fo [Hz] 

FIG. 6: The maximal amplitude of the GWs from Q-ball formation at the frequency ranging from 
1 Hz to 10^ Hz in each case is shown [Case A (red solid line), Case B (green dashed line) and 
Case C (blue dotted line)]. The horizontal axis represents the frequency, /o, and the vertical axis 
represents the present density parameter of the GWs, r^Q^*^^. The parameters are taken to be 
r/ ~ l,a ~ 1, and /i ~ ~ \K\ ~ 0.1, and the effective relativistic degrees of freedom of MSSM 
are assumed. 

production from Q-ball decay, the next-to-lightest supersymmetric particle (NLSP) decay 
may spoil the success of the BBN in some cases since the hadronic decay product of NLSP 
would destroy the light elements [3l| . 

Here we comment on the present baryon and lepton asymmetry. In the situation where the 
GWs from the Q-ball formation might be detected, both present baryon/lepton asymmetry 
and radiation are generated by the Q-ball decay. The produced baryon/lepton asymmetry 
from the Q-ball decay would be rather large and be of the order of unity unless the parameter 
in the A-term, a^, is strongly suppressed j32|. In the case of Q-ball with baryonic charge, it 
is far beyond the experimental bound on the present baryon asymmetry. Even in the case of 
L-balls, that is, Q-balls with lepton charges but without baryon charge, the situation does 
not change since the decay temperature of Q-ball is about 500 GeV and hence the lepton 
asymmetry is converted to baryon asymmetry by sphaleron process [s^ that conserves B — L 
charge. The way to avoid such large baryon/lepton asymmetry in this scenario is that the 
Q-balls are made of the AD field with B — L = 0. In this case, however, we need other 
baryogenesis mechanisms. 
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IV. DISCUSSION AND CONCLUSIONS 



In this paper, we have discussed the detectabihty of the GWs from the Q-ball formation. 
At the Q-ball formation, Q-balls with large Q can produce a large amount of GWs. However, 
such Q-balls decay slowly and they may dominate the energy density of the Universe so that 
GWs are significantly diluted. Therefore the detectabihty of the GWs is determined by 
these two competing effects. 

We have concentrated on the case of thermal log type Q-balls. Such Q-balls do not 
dominate the energy density of the Universe until the dominant potential for the AD field 
changes, after that their properties are changed. We have shown that in the gauge-mediated 
SUSY-breaking model, if the reheating temperature is Tr ^ 10^° GeV and the initial field 
value of the AD field is 0osc — Mq with 7713/2 — lOGeV and Mp ^ 10^ GeV, the present 
density parameter of the GWs from the Q-ball formation can be as large as f^^w — IQ-^^ 
and their frequency is /o — 10 Hz. Thus, it is difficult but not impossible to detect them 
by next-generation gravitational detectors like DECIGO or BBO, but the parameter region 
for detectable GWs is very small. Moreover, we have shown that it is almost impossible 
to detect GWs from Q-ball formation in other cases. In other cases when the thermal 
logarithmic potential drives the Q-ball formation, though the present amounts of the GWs 
from the Q-ball formation can be as large as f^cw — lO"®, the frequencies of such GWs are 
turned out to be very high. In the cases where zero-temperature potential terms drive the 
Q-ball formation, the present amount of the GWs from the Q-ball formation is very small 
due to the large dilution. Thus, the identification of such GWs may determine the decay 
rate of inflaton or the initial condition of the AD mechanism. In Table IIIIl we show the 
minimal frequency of the GWs that satisfy f^cw ^ 10-^^ in each thermal dominated case, 
and the maximal amplitudes of GWs in each zero-temperature case are given in Table IIVI 



Dominant term in the potential 


The smallest frequency /o that satisfies ^^q^v > 


Vihermal ^ VgraMiK > 0) 


10^-4 Hz 


^thermal Vgra.v{K < 0) 


10^-3 Hz 


^thermal ^ Vgra.v{K > 0) =^ V^auge 


10^-2 Hz 


thermal ^ ^augc 


10^-3 Hz 



TABLE IIL The minimal frequency with ^Q-y^ > 10 for thermal log type Q'ba-lls. 



Dominant term in the potential 


The maximal density parameter ^qy/ 


^grav(i^ < 0) 


10-25 


^augc 


10-21 


^gra,v{K > 0) =Y- V^augc 


10-24 


^rav(-^ > 0) =^ thermal ^ ^augc 


10-24 


thermal (c < 0) =^ V^augc 


10-24 



TABLE IV: The maximal density parameter for Q-balls with the zero-temperature potentials. 

We would like to comment on the difficulty in realizing the successful parameter region. 
One difficulty is the NLSP decay that may spoil the successful BBN when 1713/2 — 10 
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GeV, since the hadronic energy release from NLSP would cause dissociation process of light 
elements (sij. However, it can be avoided, for example, if the mass of NLSP is heavy enough 
to decay quickly. Another is the baryogenesis. One may wonder whether the present baryon 
asymmetry can be explained simultaneously for such parameter region for the detection of 
the GWs. Generally speaking, including the present case, the amount of produced baryon 
asymmetry is typically large for the case that AD condensates or Q-balls (almost) dominate 
the energy density of the Universe so that the present radiations and baryons are attributed 
to their decays. This is simply because the number densities of radiations and baryons are of 
the same order unless the (CP- violating) A-terms are suppressed by some symmetry. Thus, 
once the GWs from the Q-ball formation are detected, we have the following two possibilities. 
In the case that such Q-balls are responsible for the present baryon asymmetry, the A-terms 
are suppressed by symmetry reason. The second option is that Q-balls are irrelevant for 
baryogenesis, which is realized for the AD fields with B — L = 0. 

Although the detailed numerical calculations are required to compute the spectrum of 
GWs associated with Q-ball formation, such GWs may be differentiated from the GWs 
generated during inflation by their spectrum and from the astrophysical origin like POP III 



stars by a non-Gaussianity test [3J]. On the other hand, a first order phase transition in the 
early Universe would produce similar spectrum of GWs. However, in our case, the gravitino 
mass must be around 10 GeV for the detection of the GWs from the Q-ball formation. As 
mentioned above, such a gravitino mass may induce the NLSP decay problem. Therefore, 
if collider experiments could determine the gravitino mass by measuring the lifetime of the 



NLSP [35|, that would provide complemental information or even rule out this scenario. 

Although the identification is rather difficult, it is true that the detection of such GWs 
by DECIGO or BBO gives us information of the early Universe and the physics in the high 
energy scale, for example, it may suggest that the AD mechanism with B — L = fiat 
directions is favored.^ 

Note added in proof: Recently Ref. j4l| claims that Q balls may survive even after 
Vgrstv with K > dominates the potential. It is true that a (thin- wall) Q-ball solution exists 
for such a potential, but it is still unclear whether such a configuration is realized in the 
expanding Universe. Furthermore, even if this is the case, the dilution factor may become 
larger so that the detectability of the GWs from Q-ball formation gets even worse. Thus, 
our conclusions are unchanged. 
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Appendix A: charge evaporation from Q-balls 



As mentioned in Sec. Ill C\ Q-balls can release their charges via the evaporation and the 
diffusion effects in the thermal bath. In this Appendix, we give the condition that Q-balls 
are evaporated out by these effects before their decays. 

At finite temperature, Q-balls are surrounded by charged free particles. The minimum 
of free energy is achieved when all charges are distributed in the form of free particles in the 
thermal plasma, which induces the charge evaporation from the surfaces of Q-balls. 

Charge emission of Q-ball consists of the two processes, evaporation at low temperature 
and diffusion at high temperature. When the difference between the chemical potential of 
the plasma (fip) and that of Q-ball (/iq) is small, chemical equilibrium is almost achieved 
and hence the charge evaporation from the Q-ball surfaces is small. Instead, the charges in 
the plasmas around the Q-balls are taken away via the diffusion process. The diffusion rate 
is given by 

Tdiff = -47radifrT. (Al) 



Here a numerical factor a^is is estimated as = 4 — 6 for quarks and squarks [37|, l38| , 



and adiff = 100 — 380 for leptons and sleptons [39 . 

On the other hand, when the difference between the chemical potential of the plasma and 
that of Q-ball is large, fip <^ fig, the evaporation process from the Q-ball surface is active. 
The evaporation rate is estimated as 

r_p = ^ ^ -4vri?Je(/iQ - (A2) 

where 

1 for T > m^, 

tV , (A3) 
— for i < rrij). 

Then, depending on the type of Q-ball, we obtain the evaporation rate, 

T 

— 2^/271 ^Q^^^—j^ for the thermal log type, 

—SttC, j-777 for the gravity — mediated type, (A4) 

m^i3/2)\K\ 

— 2\/27r^Q^/^-— — for the gauge — mediated type. 

Mp 

Here we have used Eqs. (!34l) . (14T!) and (143|) and the fact hq ~ u. The rate of the charge 
transportation is determined by mindFevapl, iTdiff |}, because both the processes are necessary 
to strip the charges from the Q-balls. We define the transition temperature as Ttr, at which 
the diffusion process is more effective than the evaporation process. Hereafter we examine 
the condition for each type that the Q-balls are evaporated out before their decay. 



r 

evap 



1. Thermal log type 

We consider the case where thermal dominates the potential. Since the Q-ball transforma- 
tion temperature satisfies > in all the cases and hence Fevap > Tdiff is always satisfied. 
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we have only to consider the diffusion process. When the condition H = — iFdifrl is satisfied 

at T > Tc, Q-balls are evaporated out before their transformation, which is realized if the 
following condition is satisfied: 

47radifr T|Mg 



Q < 



47radiff Mg 



for Tn < Tc, 
for Tr>Tc. 



Here we have used the temperature dependence of the Hubble parameter: 

A 



for T>T, 



H 



4/^ TIMg 

yli/2_ for Tr<T<T, 



Mg 



R- 



(A5) 



(A6) 



2. Gravity-mediated type 



Next we consider the case where V^rav(2) dominates the potential. In this case, the tran- 
sition temperature Ttr is given by 

Ttv ^ rn^{3/2)- (A7) 

Here we have approximated l-fTladig ~ 1. 

In the case where there is Q-ball dominated era, the Hubble parameter decreases with 
the temperature as 

^ ^ for T>Tr, 
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All' TIMg 



^1/2 



for Tdom <T <Ti 



Mg ^ 

^1/2 dom T<T<T, 



R, 



M, 
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T 



for Tdcc <T <T„ 



(A8) 



A'/^^ for T<Tdec. 



M, 



G 



Here Tp = {T^^^TdomY^^ ■ Then min{|revap|, |rdiff|}/(<5-ff) becomes the highest at T ~ 
'^<^(3/2)- Thus, the Q-balls are evaporated out before their decays if the following condi- 
tion is satisfied: 



Q < < 



47radiff4' T^Mg 

^0(3/2) ^1^3/2) 

47radiff Mg 



diff 



M, 
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41/2 1/2 ^1/2 

^(/>(3/2) ''>(3/2)^dom 



for Tr < m<^(3/2), 

for Tdom < "10(3/2) < Tr, 
for m^(3/2) < Tdom- 



(A9) 



Here the subscript "0(3/2)" indicates that A is estimated at T = m(^(3/2). 
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3. Gauge-mediated case 



Here we consider the case where Vgauge dominates the potential. In this case, the transition 
temperature T^j. is given by 



tr 



1 (v^adiffm2M^Q-i/4)i/3 fo^ T,, < 



for Ttr > m^, 



(AlO) 



These two temperatures coincide for Q ^ Qcr = a^{Mp/m^Y. 

Using the similar argument above, we have the condition that the Q-balls are evaporated 
out before their decays. 



f / ^ V 4/3 / „ X 4/3 
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(All) 
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TLMg] 
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for Tdcc < ^tr < Tp 



In conclusion, Q-balls cannot be evaporated out before their decays unless the charge Q 
inside a Q-ball is small enough. 

Appendix B: The case with the zero-temperature potential 

In this appendix, we show that the GWs from the Q-ball formation with zero-temperature 
potential are too small to be detected. 



1. Gravity- mediated type 

First we consider the gravity-mediated or the "new" type Q-ball. This type is real- 
ized for the gravity or anomaly mediated SUSY breaking model with < or for the 
gauge-mediated SUSY-breaking model with V^rav(2)(-^' < 0). In this type, the potential is 
dominated by 

i + ^log(i-L 



^grav(2) — ''^0(3/2) 



M, 



G 



l$l 



(Bl) 



and from Eqs. 
estimates as 



and fl25p . the parameters associated with the Q-ball formation are 



n '-' I ""m 



K\, and ~ 



"^0(3/2) 1^1 

2a 



(B2) 
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Here, the factor a ~ 30 represents the dilution due to the cosmic expansion. Then, the 
amount of the GWs, f2Gw(^max/(7ra)), at the Q-ball formation is estimated from Eq. ( 168|1 

as ^ 

nl^{^rQm,\K\^''|{2^)) ^ ^a'^' ||. (B3) 



The baryon or lepton charge stored in a produced Q-ball is estimated as [23 



Q^^(^^] , (B4) 

where the numerical factor /3 ^ 6 x 10~^ represents the dilution due to the cosmic expansion. 
Other properties of the Q-balls are given by 

2 2 

^ \K\' a;~m<^(3/2), 

''>(3/2)l-^ I 

f\K\V^^ ~ 1 
(pQ-i^j /3^^Vosc, ~ -m^(3/2)Q. (B5) 



Then, the average density of the Q-balls can be estimated as 



V 



(B6) 



Here 77 is a numerical factor of order unity. Now we investigate the present properties of the 
GWs from the Q-ball formation. One should notice that there is no Q-ball transformation 
different from the case when the thermal logarithmic potential dominates the potential. 
Therefore, we have only to consider two cases depending on whether there is a Q-ball 
dominated era or not. Figure [7] shows the schematic time evolution of the energy density of 
each component. The Hubble parameters at the Q-ball domination (if any) and the Q-ball 
decay are given, respectively, by 



/3 ^mj] 



osc 



24710; 

Here we have used Eq. (149|) and the relations 

3r/2|ir|2M2' 

i/dom ^ ^qHr. (BIO) 

Then, we have the condition for the Q-ball domination, 

^ / 3 Y" ,f'^\K\'"' Mfmf 
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FIG. 7: The time evolution of the energy density of each component is shown schematically. While 
Prad aiid pgw decrease in proportion to a~^, pinf in the inflaton osciUation era, and pq decrease in 
proportion to 

The present density parameter of the GWs from the Q-ball formation in each case is given 
by 



^^GW — 



192 Vvr; 7r3 Mt!^ «o 



— ) - — ' ^ — '- TTT— (with the Q-bah domination). 



12 M}."^mT «o 



(B12) 

osc-'R "eq (^^jthout the Q-bah domiuatiou) . 



G 



Here we have used Eqs. fl83|) . (lB2p . flB7l) and flB8l) . Then, from Eq. fIBlip . we have the upper 
bound on the present density parameter of GWs, 

- 192 yn) TT^ M'f ao ^ ^ 

which is too smaU to be detected by the future detectors. 

Here we comment on the possibihty of the evaporation of Q-balls before their decay. From 
Eq. ( 1A9I) . we find that the Q-balls are evaporated out before their decays if the following 
condition is satisfied: 

^ 47ra Mg 

Q < -Y7^ B14 



>(3/2) 



which yields, using Eq. flB4p . 



Ana 



1/2 



<(3/2)^G ' ^ lO^^GeV. (B15) 



.1/2 n j ">(3/2) 
. </.(3/2)/^ 

Such a small initial field value does not produce significant amounts of GWs and have the 
present amplitude of GWs is rather small. Thus, our conclusion is unchanged even if we 
take into account of the possibility that the Q-balls are evaporated out before their decays. 
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2. Gauge-mediated type 

Next we consider the gauge-mediated type Q-balls. In this type, the potential is domi 
nated by 



Vi.„,e=iMJ(^logi^l . (B16) 

This potential is quite similar to the thermal logarithmic potential. Then, the parameters 
associated with the Q-ball formation are given by 

*=L._,3M| ML (BIT) 



a 



which yield the amount of the GWs fiGw(^max/(7ra)) at the Q-ball formation from Eq. (l68ll . 



G 

1/2 

Other properties of the produced Q-balls are obtained by replacing ag Tqsc by Mp in Eq. 



m 

Now we evaluate the present properties of the GWs from the Q-ball formation. Note 
that there is no Q-ball transformation in this type and hence the cosmic history depends on 
whether their is Q-ball dominated era or not, as is the case with the gravity-mediated type 
(Appendix IB ip . The Hubble parameters at the Q-ball domination (if any) and the Q-ball 
decay are, respectively, given by 

i^dec - ^(e/3) (B20) 



osc 



which give the condition for the Q-ball domination 

1/4 



Here we have used Eqs. (H9l) and (IBIOI) . Using Eq. ( I83l) . the present density parameter of 
GWs from Q-ball formation is estimated as 

[^\ (6/3)-5/6aV^' ./f./, — (with the Q-ball domination), 

_ I 216 VV2y M^f(t)%l «o 

"GW — \ ^8/3 ^ ^ 014/32^4/3 ^ 

— —{ej3)A^J(^ — °7/3 ^4/3 ~ (without the Q-ball domination). 

(B22) 

The condition that (pg = (e/3)^/^0osc > Ms = Mp/m^ yields the upper bound on I^q^' 

1 /q 2\2/3 4/3 
f]0,max ^^(^\ (^^)-5/6^4 4/3^^ ^ ^q-21_ ^333) 

«w 216 V 72/ Mj/'ao ^ ' 
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Again, the GWs generated by the formation of the gauge-mediated type Q-balls are too 
small to be detected by the next-generation GW detectors. 

Finally, we give the condition that the Q-balls are evaporated out before their decays, 

fM \^/^ 

Q < [j^) ^ 0OSC < MfM-"\ (B24) 

which yields 0osc ^ 10^ GeV since Mp < 10^° GeV. Such a small initial field value does 
not produce significant amounts of GWs and hence the present amplitude of GWs is rather 
small. Thus, our conclusions are unchanged even if we take the Q-ball evaporation into 
account. 



3. Delayed type 

Next, we consider the delayed type Q-balls. In the gauge- mediated SUSY-breaking model, 
the effective potential for the AD fields is given by the summation of V'grav2, ^augc, and 
thermal- If the AD field starts its oscillation where the potential is dominated by V'grav2 with 
-ff > 0, Q-balls are not formed and the AD field falls down along the potential. Then, thermal 
or Vgauge dominates the potential of the AD field at a critical temperature or a critical field 
value, which induces the formation of Q-balls, called the delayed type Q-balls. Thus, there 
are two cases depending which potential term is responsible for the Q-ball formation. 



0- ^augc driven Q -hall formation 

First we consider the case where Vgauge drives the Q-ball formation. In this case, when 
= M|,/m3/2, Q-balls are form 

given by 



eg = M|./m3/2, Q-balls are formed. The Hubble parameter at the Q-ball formation is 



^ (B25) 

t'osc 

where 0osc is the field value at the onset of the oscillation of the AD field and we have 
used the relation (j) (x H during the oscillation of the AD field. The condition that Vgauge 
dominates the potential before thermal is given by Tr < Mp. From Eqs. and ( j25l) . the 
parameters associated with the Q-ball formation are given by 



/5g,.^^ and ^^mj/,, (B26) 



which yield the amount of the GWs, f2Qw(^max/ (^o)); at the Q-ball formation from Eq. fl68|l . 
The baryon or lepton charges stored in a Q-ball are given by [Isf 
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Here the numerical factor /3 [defined in (133|) ] is almost unity since the Q-balls are formed very 
quickly so that the cosmic expansion is negligible. Other properties of Q-balls are estimated 
as, 

^ M^gV4 ^ i^L^ ^ l^^M^g^/^ (B29) 

The energy density of Q-balls at the Q-ball formation is given by 

p*Q^EQ{klJa^)^M^p. (B30) 

Then, we evaluate the present properties of the GWs from the Q-ball formation. Again, 
there is no Q-ball transformation and hence the cosmic history depends on whether there is 
the Q-ball dominated era or not. The Hubble parameters at the Q-ball domination (if any) 
and the Q-ball decay are, respectively, given by 

4 1/2 ,4 ry2 

H,or. = ' (B31) 

m|/2 



Here we have used the relations, 



ffdom ^ (I'^Hr. (B34) 

Thus, the condition for the Q-ball domination is given by 

From Eq. f l83p . the present density parameter of GWs from Q-ball formation is given by 



216 V 72/ M'/' ao 



(with the Q-ball domination). 



A^^&^M^Oeq ^^^^^ 

54 Mi^/^m^/, ao 



(without the Q-ball domination). 



Considering the condition, that 1713/2 < lOGeV for the gauge- mediated SUSY-breaking 



model, we have the constraint on f^GW; 



fi°,w < 10-'^ (B37) 

which is too small for the detection. This conclusion applies also for the second Q-ball 
formation discussed in Sec. HHP 31 
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Here we comment on the Q-ball evaporation. From Eqs. (lAlip . the Q-balls are evapo- 
rated out before their decays if the following condition is satisfied: 

Q < {^)'"' ^ M, < M^!'ml',t. (B38) 
Thus, we have the inequality with respect to the amount of the GWs at the Q-ball formation, 

Note that this amount of GWs is further diluted at least after matter-radiation equality. 
Therefore we conclude that even if we take into account the possibility that the Q-balls 
are evaporated out before their decays, the amount of the present GWs is too small to be 
detected by the next-generation detectors. 



b- ^thermal driven Q -ball formation 

Next we consider the case where thermal drives the Q-ball formation. In this case, when 
T = and 0* ^ agT^ /'m3/2, Q-balls are formed. Note that the reheating takes place after 
the Q-ball formation because we have the following inequality: 



In this case, we have two constraints on the reheating temperature. One is the condition 
that V'grav2 dominates Vthermai at the onset of the AD field oscillation, 

Tr < ^ ^0osc = T^d''- (B41) 

Another condition is that Vthermai dominates Vgauge at the Q-ball formation, 

Tn > ^^Mj. ^ T^. (B42) 
The properties associated with the Q-ball formation are given by 



^r^^ and ^^ml/„ (B43) 



which gives the amount of the GWs, f2Q^(A;max/(7ra)), at the Q-ball formation. 
Here we have used the following relations: 



- " # (07^ ' ^^^^^ 
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The baryon or lepton charges stored in a Q-ball are given by |23l 



^ 4 / X 2 / ^ X 4 



\a/ TJ \(poscJ V^3/2/ 

Here the numerical factor /3 is almost unity since the Q-balls are formed very quickly so that 
the cosmic expansion is negligible at the formation in this case. Other properties of Q-balls 
are evaluated as 

n ~ UJ ~ 



^ ay2TgV^ Eq ^ ^^al/'TQ-'/\ (B47) 



At the Q-ball formation, the energy density of Q-balls is given by 

and it decreases with the temperature as 

PQ oc A^T^ (B49) 

Next, we follow the cosmic history after the Q-ball formation and evaluate the present 
properties of the GWs from the Q-ball formation in this case. The properties of Q-balls 
are changed into those of the gauge-mediated type Q-balls, when ^auge dominates thermal 
at T ^ T(7 = ag^^'^Mp. In the same way as the cases discussed in Sec. IIII Ct and Sec. 
nil D[ where the thermal logarithmic potential dominates the effective potential, we have 
four possibilities of the cosmic history after the Q-ball formation. They are characterized 
by the critical temperatures, and T^. In this case, they are summarized in Table IVl 



Case 1 


Tn > max{Tf (0osc),T^'(</'osc)} 


Case 2 


Ti{<Posc) <Tr< T^^{<I)osc) 


Case 3 


Ti{(t>osc) < Tr < T^H<l>osc) 


Case 4 


Tr < min{r^i(0osc),71^(</'osc)} 



TABLE V: The conditions of four cases of the cosmic history 

The critical temperatures are given by 

T'r' = a-^'^Mp, (B50) 



Here we have used the Hubble parameters at the Q-ball domination (if any) and the Q-ball 
decay given by 
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Then, from Eq. (183|) the present density parameter of GWs from Q-ball formation is esti- 
mated as 



0° ~ 
"GW — 



1 /Q 2\ 2/3 2/3.1/6 T^/^TO^/^ 



54 Mi°/^m^/,^ 



(B54) 



(Case 2,4). 



The condition Eq. (]B4ip yields the upper hmit of the present amount of the GWs in Case 
1 and Case 3, 

which is too small to be detected. Here we have used the fact that Ug ~ O.l,0osc < 
Mg, 1713/2 < 10GeV,MF > lO'^GeV, and tteg/ao ~ 3 x IQ-^ and g^,s are evaluated in 
the context of MSSM. Moreover, in Case 2 and Case 4, the constraint on the reheating 
temperature Tr < strongly constrains the upper bound of ^%iy and hence the present 
amount of the GWs cannot be larger than that in Case 1 and Case 3. Therefore we conclude 
that the GWs generated by the formation of this type of Q-balls are too small to be detected 
by the next-generation GW detectors. 

Here we comment on the effect of the Q-ball evaporation. From Eq. f lAlip . the Q-balls 



are evaporated out before their decays if the following condition is satisfied: 

\MfJ ag M^^^Mp'^ 

which gives the inequality for the amount of the GWs at the Q-ball formation, 

1 TTl^ d)'^ 

Considering the fact that 1713/2 ^10 GeV, 0osc ^ Mq, and Mp > 10'' GeV, we have the 
upper limit on fiQy/, 

< 10-^^ (B58) 

which is further diluted at least after matter-radiation equality. Therefore we conclude that 
even if we take into account the possibility that the Q-balls are evaporated out before their 
decays, the amount of the present GWs is too small to be detected by the next-generation 
detectors. 



4. Negative thermal log type 

Finally, we consider the case with negative thermal logarithmic potential. So far, we 
assumed that the contribution of the thermal logarithmic potential is positive. However, it 
is possible for the term to be negative and hence there can be another type of Q-balls [40|. 

We then consider this type of Q-balls. If the temperature after inflation is sufficiently 
high, the AD field is trapped in the potential minimum, ~ [ag 

rp2j^n-3y/(n-i)^ determined 
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by the balance between the nonrenormahzable F-term and the negative thermal log term, 
rather than the negative Hubble mass term. In the gravity or anomaly mediated SUSY- 
breaking model, this potential minimum exists until the thermal correction to the potential 
turns to the thermal mass term. At that time, the AD field starts oscillating around the 
origin by the thermal mass term so that it decays quickly. Thus, Q-balls are not formed in 
this case. On the other hand, in the gauge-mediated SUSY-breaking model, this potential 

Mp- In this case, the AD field starts oscillating from 



-1/2 



so that Q-balls are formed. Hereafter, we consider such 



minimum vanishes when T = ag 
0OSC ^ (M2T2M"-3)V(n-i) by V„ 
a type of Q-balls. 

In this case, the parameters associated with the Q-ball formation are given by. 



'■'max 



3M) 



2 0L' 



/3, 



^ and 



for Tf? > a, 



-1/2 



The properties of the Q-balls are estimated as 



Mi 



^24/2 MgTI 



for Tr < a^^^^Mp. 



(B59) 



Q 



^ 0OSC, and 



(B60) 



Then, the amounts of the GWs, ^gW' Q-ball formation are given by 

.,2 / i \ 2 



^GW — 



54 A \Mg 



[ 54^2 \Mg) Xat'^Mp 



for Tr > ag ^^"^Mpi 



for Tr < a^^^^Mp. 



(B61) 



Note that since the Q-balls are formed quickly, there are no dilution factors. 

Next, we follow the cosmic history after the Q-ball formation and evaluate the present 
properties of the GWs from the Q-ball formation. In this case, there are four possibilities 
of the cosmic history. They are classified by two criterions. One is whether the Q-ball 
dominated era exists or not and the other depends which first takes place, the reheating 
after the infiaton decay or the beginning of the oscillation of the AD field. 

The Hubble parameters at the Q-ball domination and at the Q-ball decay, H^^c are, 
respectively, given by 



H. 



dec 



9A'J' Mg 

^843/2 J^Q 

^ 9AI MgM^ 
48 6^ ' 

^ Vosc 



for Tr > a'g^'^Mp 
for Tr < a-g^'^Mp 



(B62) 



(B63) 
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which gives the condition for the Q-ball domination, 



> < 




3/5 ^^^G ^^^F 



8/5 ^3/10 



for Tr > a\ 



for Tr < a. 



-1/2 



(B64) 



-1/2 



The present density parameter of the GWs from the Q-ball formation is given by 



oo ^ 
"gw — 



216 

„2 



3r 



2/3 



8/3 



"eq 



72/ M 



4/3 ,4/3 



G 



with the Q-ball domination. 



54 A V^^G 



ao 



— 1/2 

without the Q-ball domination and Tji > ag M 



54^2 \M, 



R 



20/3 



"'eg 

ao 



without the Q-ball domination and Tr < a. 



-1/2 



(B65) 



Then, the upper bound on ^Q-i^ can be obtained by the same consideration as the case 
discussed in Appendix IB 21 and hence are given by ^q^**^ ~ 10"^'^. Moreover, the possibihty 
of the Q-ball evaporation does not change the conclusion as is the case in Appendix IB 2[ 
Therefore we conclude that the GWs generated by the formation of this type of Q-balls are 
too small to be detected by the next-generation GW detectors. 

We conclude that the GWs from the Q-ball formation can hardly be detected by the future 
detectors in the case where the oscillation of the AD field is driven by the zero-temperature 
potential terms. 
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